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ABSTRACT 

Oxygen abundances of 67 dwarf stars in the metallicity range —1.6 < [Fe/H] < —0.4 are derived from 
a non-LTE analysis of the 777 nm Oi triplet lines. These stars have precise atmospheric parameters 
measured by Nissen and Schuster, who find that they separate into three groups based on their 
kinematics and a-element (Mg, Si, Ca, Ti) abundances: thick-disk, high-a halo, and low-a halo. We 
find the oxygen abundance trends of thick-disk and high-a halo stars very similar. The low-a stars 
show a larger star-to-star scatter in [O/Fe] at a given [Fe/H] and have systematically lower oxygen 
abundances compared to the other two groups. Thus, we find the behavior of oxygen abundances in 
these groups of stars similar to that of the a elements. We use previously published oxygen abundance 
data of disk and very metal-poor halo stars to present an overall view (—2.3 < [Fe/H] < +0.3) of 
oxygen abundance trends of stars in the solar neighborhood. Two field halo dwarf stars stand out in 
their O and Na abundances. Both G53-41 and G150-40 have very low oxygen and very high sodium 
abundances, which are key signatures of the abundance anomalies observed in globular cluster (GC) 
stars. Therefore, they are likely field halo stars born in GCs. If true, we estimate that at least 3 ± 2 % 
of the local field metal-poor star population was born in GCs. 

Subject headings: stars: abundances — stars: Population II — Galaxy: disk — Galaxy: halo 



1. INTRODUCTION 

Important clues to understand the formation and evo- 
lution of the Milky Way's halo and disk components, 
as well as any possible connections between them, are 
imprinted in the photospheric chemical composition of 
FGK-type dwarf stars. These objects have probably re- 
tained the chemical composition of the gas from which 
they formed, thus being excellent tracers of Galactic 
chemical evolution (GCE). Moreover, their long lifetimes, 
in particular for the G and K types, allow us to probe 
the Milky Way's GCE over many billions of years. Con- 
sequently, combined with kinematics and information on 
stellar ages, chemical composition analyses of FGK dwarf 
stars have the potential to be useful for reconstructing 
the history of our Galaxy. 

The simplest picture for the Galactic halo formation 
involves a monolithic collapse which leads to a halo star 
population showing a strong correlation bet ween orbital 
eccen tricity an overall metal abundance (|Eggen et al.l 
1962). Although historically important, this model has 
long been known to be incomplete. Metallicity determi- 
nat ions of giant s t ars in globular clusters, for example, 
led iSearle fc Zinnl (|1978f ) to conclude that the halo was 
formed in a more "chaotic" fashion. Indeed, state-of-the- 
art simulations show that the properties of the stellar 
component of galactic halos could be heavily influenced 
by m erging events that occur as the galaxy assembles 
(e.g.. lAbadi et al.l[200l IGuo fc White! 120081: [Read et al.1 
120081 : iStewart et al.l l2008t iScannapieco et al.ll2009h . In 
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these scenarios, rather than consisting of a single sim- 
ple evolving population, the halo is expected to contain 
sub-structures as remnants of its formation history. The 
discovery and detailed characterization of hal o streams 
and tidal debris h e avily support this idea fe-g.-lMaiewskil 
119931 : iHelmil 120081 : iKlementlboiOt iMaiewski et aLl l20121 . 

Very strong evidence for halo sub-structures in the so- 
lar neighborhood has been recently provided in a se- 
ries of papers by INissen fc Schusterl pcM |20H . and 
iSchuster et al.l (|2012ft . From a detailed spectroscopic 
analysis of 94 dwarf stars in the [Fe /H] range from 
-1.6 to -0.4 fflNissen fc Schusterl lj2O10L hereafter NS10) 
found that stars with halo kinematics separate into two 
groups based on their a-element abundances (in their 
case quantified by the average abundance of Mg, Si, Ca, 
and Ti). NS10 argue that the halo stars in the low-a 
group could have been accreted from satellite galaxies, 
possibly u) Cen. Their abundance analysis of heavier ele- 
ments, particularly Na and Ba/Y, however, showed that 
the <jj Cen and low-a halo star connection is weak, un- 
less chemical evolution within the satellite galaxy was 
different for its inner and outer regions, an idea sup- 
ported by the obser vation of an abundance gradient in 
a dwarf galaxy (cf. INissen fc Schusterl l2?)ll . Finally, 
ISchuster et al.l (|2012f ) show that the low-a halo stars are 
about 2-3 Gyr younger than the high-a halo stars and 
that these two groups exhibit different orbital properties, 
with the low-a stars having very eccentric orbits, larger 
r max (maximum distance from the Galactic center), and 
larger z max (maximum distance from the Galactic disk). 

A very important chemical element missing from the 
NS10 paper series is oxygen. As the third most abundant 
element in the universe and in stellar atmospheres, after 

4 In this work we use the standard definitions: [X/Y] = 
log(JVx/iVV) - log^x/AV)©, and Ax = log(Nx/Nn) + 12, where 
Nx is the number density of element X. 
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H and He, and having one of the best identified produc- 
tion sites of all elements as well as reliable supernovae 
yields, oxygen is crucial for GCE studies. Furthermore, 
oxygen is a key element in the investigation of abun- 
dance variations in globular clusters (GCs). Stars with 
enhanced Na are known to be depleted in O, i.e., they 
follow the we ll-known oxygen-sod i um anti-correlation 
in G Cs (e.g.. IGratton et all l200i fCohen fe Melended 
[2005t lYong et al.H2005HAlves-Brito et al.ll2012ri . A' num- 
ber of recent studies have investigated the contribu- 
tion of GCs to the build-up of the field halo population 
(lYong et al.ll2008HCarretta et alJl20intlMartell fc Grebel 
l2010HMartell et al.ll2011l ). but none of them have found 
field stars with both high Na and low O abundances. 
Since NS10 have already studied Na in their sample of 
metal-poor stars, the addition of oxygen allows us to as- 
sess to which level the halo field has been contaminated 
by stars formed in GCs. 

Determining reliable oxygen abundances in metal-poor 
dwarf stars is not a straightforward task. Few spectral 
features due to oxygen are available in the visible spec- 
trum and all are affected by a number of model uncertain- 
ties or severe line blending. Our past experience success- 
fully deriving oxygen abundances from a restricted non- 
LTE analysis of the 777 nm O I triplet ( Ramirez et al. 
[20071 hereafter R07) now allows us to infer them. In this 
work, we derive oxygen abundances for as many as pos- 
sible of the stars in the NS10 study in order to better 
understand the nature of the two distinct halo popula- 
tions in the solar neighborhood. 

2. SAMPLE AND SPECTROSCOPIC DATA 

Given the careful sample selection and high precision 
of the stellar parameters and elemental abundances de- 
rived by NS10, we adopted their sample in our work, as 
stated above. We collected high-quality spectra of the 
777 nm region for as many of these stars as possible. We 
started by searching for data in our own spectral libraries 
and then in publicly available data archives. We comple- 
mented this data set with new observations, as described 
below. 

First, we used data from the R07 work on oxygen abun- 
dances in nearby stars. Nearly all of these spectra were 
acquired with the R. G. Tull coude spectrograph on the 
2.7 m Telescope at McDonald Observatory, and reduced 
in the standard manner using IRAF0 as described in 
R07. One spectrum from the High Resolution Spec- 
trograph (HRS) on the Hobby-Eberly Telescope (HET), 
reduced also as in R07, and one spectrum from the 
VLT-UVES POP (P aranal Observatory Project) library 
(jBagnulo et al.ll2003[ ) were also used. The McD-Tull data 
have a spectral resolution R = A/AA ~ 60, 000 while the 
HET-HRS and VLT-UVES spectra have R ~ 120 000 
and R ~ 80 000, respectively. 

Then, we se arched for spectr a taken with the HIRES 
spectrograph ()Vogt et al.lll994D at the Keck I Telescope. 
We found 11 stars available in the Keck Observatory 
data archive covering the 777 nm O I triplet. These spec- 
tra were reduced using MAKEE0 a data reduction tool 

5 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation — http://iraf.noao.edu 

6 http:/ /spider. ipac.caltech.edu/staff/tab/makee 



developed by T. A. Barlow specifically for reduction of 
Keck-HIRES spectra. In some cases, we re-reduced the 
archive spectra with MAKEE by fine-tuning the extrac- 
tion parameters in order to improve the results. Most of 
the HIRES spectra have R ~ 67 000, but some of them 
were taken at R ~ 50 000 or R ~ 100 000. 

We found 27 of the NS10 stars in the R07 work and 
11 more in the Keck-HIRES archive. In order to build 
a more statistically significant sample, we performed 
spectroscopic observations of 24 additional stars using 
the MIKE spectrograph on the 6.5 m Magellan/Clay 
Telescope at Las Campanas Observatory. We acquired 
the data for these stars in four observing runs in July, 
September, and November of 2011, as well as in February 
of 2012. We used the MIKE standard setting with the 
narrowest slit (0.35arcsec width), which provides com- 
plete wavelength coverage from 3350 to 9500 A, includ- 
ing the 777 nm oxygen triplet, at a sp ectral resolution 
R = 65 000 in the oxygen triplet region (Bernstei n et al.l 
2003). The signal-to-noise ratios of these spectra vary 
from star to star from about 100 to 500 at 777 nm. The 
MIKE spectra were reduced using the CarnegiePython 
pipeline^ One more star (G150-40) was observed at 
McDonald Observatory in April of 2012. Its spectrum 
(R ~ 60 000, S/N ~ 150) was reduced as in R07. 

The high quality of our spectroscopic data allows a very 
precise measurement of the equivalent widths (EWs) of 
the three lines of the oxygen triplet. Based on the spec- 
tral resolution and typical signal-to-noise ratio of most 
of our data (R ~ 60 000, S/N ~ 200), we estimate an 
err or of about .7 mA for each EW measurement (us- 
ing [CiyTej|[l988| formula). The actual EW error is likely 
larger due to uncertainties in the continuum determina- 
tion, but tests made measuring the EW of a number 
of lines while varying slightly the continuum location 
showed that these errors are unlikely to be larger than 
about 1.0 mA. Thus, we estimate our EW errors to be 
about 1.2 mA. 

Three of the stars from the NS10 work, and for which 
we did not find spectra or were unable to acquire them, 
have prec ise (errors of order 2 m A) EWs of the oxygen 
triplet bv lAkerman et al.l (|2004[ ). We adopted those EW 
values in our work. Also, for G75-31, although a spec- 
trum from R07 is available, it is of relatively low quality, 
and therefore we preferred to use th e EW values for the 
star measured bv lNissen using a higher qual- 

ity VLT-UVES spectrum. For the Sun, which we use to 
transform the abundances from absolute (^o) to relative 
([O/H]), we adopted the solar EWs by R07. These values 
were obtained as the average of EWs measured in three 
solar spectra, two skylight observations and one asteroid 
reflected sunlight spectrum, which were shown to be in 
very good agreement. 

3. OXYGEN ABUNDANCES 

The EWs of the three lines of the O I triplet at 777 nm 
were measured by fitting Voigt profiles using IRAF's 
splot task, except for weak lines with relatively low 
signal-to-noise ratio spectra. The Gaussian profile is a 
very good approximation to the real line shapes consid- 
ering the spectral resolution of most of our data and the 
fact that these lines are not so strong in this metallic- 

7 http:/ /obs. carnegiescience.edu/Codc/mikc 
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ity regime. Voigt profile fits are more accurate for strong 
lines with extended wings, but they tend to confuse noise 
with wing depth for weak lines in low S/N spectra. In 
those cases a Gaussian fit often works better. Our mea- 
sured equivalent widths were then used to derive the oxy- 
gen abundances using a standard curve-of-growth (COG) 
approach, as described below. 

The abfind driver of the 2010 vers ion of the spectrum 
synthesis code MOOG (ISnedenlllQ?^ was used to com- 
pute COGs and derive the oxygen abundances from our 
EW measurements. We used the new MA RCS grid of 
model atmospheres ([Gustafsson et al.l I2008H . The stel- 
lar parameters T e g, logg, and [Fe/H] adopted are those 
derived by NS10. For a number of stars (6), NS10 pro- 
vide two sets of stellar parameters, each derived from 
a different spectrum. For all these stars, the two sets 
of parameters are in excellent agreement considering the 
observational errors. In these cases, we adopted the av- 
erages of the two sets of atmospheric parameters derived 
by NS10. 

NS10 estimate their differential errors for T c g, \ogg, 
and [Fe/H] at 30 K, 0.05 dex, and 0.03 dex, respectively. 
Using the reference stars HD 22879 and HD 76932, we can 
obtain an estimate for the oxygen abundance ( [O /H] ) er- 
ror by calculating the variation of their absolute oxygen 
abundances, Ao, for a given change in stellar parame- 
ters. This procedure is safe because these uncertainties 
are only weakly correlated. Propagated into our abun- 
dance analysis, the errors in stellar parameters translate 
into a 0.025 dex uncertainty for [O/H]. The microturbu- 
lcnt velocities, v t , were also adopted from NS10, where 
no estimate of the vt error is given. If we assume a Vt 
error of 0.2 km s -1 , which is a conservative estimate con- 
sidering the high precision of the strictly differential work 
by NS10, the uncertainty introduced to the oxygen abun- 
dance is only about 0.015 dex. The EW uncertainty of 
each of the triplet lines translates into a 0.020 dex error, 
which implies a total oxygen abundance error of about 
0.035 dex (internal error only). Model uncertainties and 
noise or instrumental defects are probably responsible for 
the larger line-to-line scatter (in the oxygen abundances 
inferred for the three of the triplet lines) seen in some of 
our sample stars. 

The spectrum synthesis code used to derive our oxy- 
gen abundances assumes local thermodynamic equilib- 
rium (LTE). However, it is well known that in the pho- 
tospheres of cool stars, both dwarfs and giants, the O I 
777 nm triplet spectral feature is formed under condi- 
tions far from the LTE approxi mation. Us i ng a well- 
justified two-level approximation, IKiselmanl (|1993l) ele- 
gantly demonstrates that the non-LTE effect is due to an 
infrared mean intensity that departs from its LTE value 
(the Planck function) at depths where the t r iplet lines 
are f ormed (see also lEriksson fc Toftl Il979t IKiselmanl 
2001). Since the radiation field is stronger in warmer 
stars and the gas densities are smaller in lower sur- 
face gravity stars, the errors due to non-LTE effects 
are larger for warmer T e g and lower logg values. An 
important [Fe/H] -dependence is also anticipated, with 
larger non-LTE effects at lower [Fe/H] due to a de- 
crease of c ontinuum opac i ty an d smaller rate of collisions 
(see, e.g., iFabbian et al.l [20091) . In the solar case, the 

8 http://www.as.utexas.edu/ chris/moog.html 



non-LTE correction to the oxygen abu ndances inferred 
from the triplet is about 0.2 dex (e.g.. IKiselmanl 119931 : 
lTakedalfl99l iRamirez et ail [20071 IFabbian et all 12009(1 . 
In warmer stars (T e g ~ 6250 K), these corrections can be 
as high as ~ 0.4 dex, which means that even in a solar 
relative analysis, differential non-LTE errors of ~ 0.2 dex 
could be introduced. Thus, non-LTE corrections to the 
LTE abundances derived with MOOG must be applied. 

A number of authors have computed non-LTE correc- 
tions for FGK-type stars of several metallicities, includ- 
ing those of our targets. There is good qualitative agree- 
ment between the various calculations available in the 
literature, but differences of order 0.1 dex between them 
are common. In this work, we use the non-LTE cor- 
rections tabulated by R07. The main reason for this 
choice is that the same corrections were applied to the 
abundances of disk stars which we will use later in Sec- 
tion 21 and similar correcti ons were applied to th e very 
metal-poor stars from the (Melend ez et al.ll2006l here- 
after M06) work, which we will also use in Section [4] To 
prevent systematic biases arising from the use of very 
different non-LTE prescriptions in our final analysis, we 
decided to make a choice that will bring the non-LTE 
oxygen abundances to approximately the same scale. 

R07 computed their non-LTE corrections using an oxy- 
gen model atom with 54 l evels and 242 transitions by 
lAllende Prieto et al.l (120031) . with a few minor improve- 
ments! IKuruczl (|1993l ) atmosphere models were em- 
ployed to determine the level populations by solving the 
rate equations with TLUSTY (|Hubenvlll988l ). Although 
this implies a certain level of inconsistency within this 
work, because the model atmospheres used to derive the 
LTE oxygen abundances are from the MARCS grid, note 
that the non-LTE corrections were deriv ed with r espect 
to LTE abundances computed also from IKuruczl (1993) 
models. It is expected that the non-LTE corrections 
are less sensitive to the choice of model atmosphere grid 
than the absolute val ues of the abundance s, either LTE 
or non- LTE. Indeed, IFabbian etaLI ([2009) have shown 
that the model-dependence in this context is only im- 
portant for stars with [Fe/H] < —2.5, i.e., stars with 
metallicities below that of our most metal-poor sample 
star. Once the non-LTE level populations were computed 
with TLUS TY, spectrum synthesi s was performed using 
SYNSPEC (IHubenv fc LanzllT99l . 

As acknowledged by R07, one of the deficiencies of their 
non-LTE calculations is the neglect of inelastic collisions 
with neutral H, which tend to bring the level popula- 
tions closer to their LTE values. R07 noticed that their 
non-LTE oxygen abundances presented systematic off- 
sets between the three lines of the triplet, and suggested 
empirically correcting for these offsets to reduce the line- 
to-line scatter. The nature of these offsets is qualitatively 
well explained by the fact that inelastic collisions with 
neutral H were ignored, but detailed calculations were 
not made to confirm this hypothesis. However, if true, 
these empirical corrections are roughly taking this effect 
into account. 

The importance of inelastic collisions with neutral H 
in non-LTE calculations is typicall y param e terize d by 
the multiplicative factor ffn to the iDrawi nl (fl968h for - 
mula, as sugg ested by iSteenbock fc Holwegerl (|1984l) . 
Interestingly, t he im proved non-LTE computations by 
IFabbian et al.l (|2009l ) . including the impact of collisions 
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with neutral H, show that their non-LTE corrections to 
the solar oxygen abundance with Sh = 1 are, albeit 
fortuitously, in good agreement with those by R07. In 
their three-dimensional hydrodynamic an alysis of oxy- 
;en line formation in the solar photosphere. IPereira et aLl 
2OO9a0) find that Sh — 1 provides an exc ellent fit to the 
observ ational data. Note, however, that Ramirez ct al. 
(2006) find that a one-dimensional static model atmo- 
sphere spectrum with Sh = 10 reproduces the 777 nm 
O I triplet line profile better than one with Sn = 1 in 
the case of the moderately metal-poor star BD +174708, 
although the derived oxygen abundance in that case ap- 
pears too high. 

We measured oxygen abundances, Ao, for each of the 
three triplet lines. The solar values were used to con- 
vert these quantities into [O/H] by averaging the three 
relative abundances. The same solar-differential line-by- 
line procedure was employed to derive the LTE and non- 
LTE relative abundances. Finally, [O/Fe] values were 
determined using the very precise (~ 0.03 dex error) iron 
abundances inferred by NS10: [0/Fe] = [0/H]-[Fe/H]. 
The mean error in the non-LTE [O/Fe] abundance ra- 
tios derived in this work is 0.05 dex. 

Our derived oxygen abundances, both in LTE and non- 
LTE, as well as the stellar parameters by NS10, their 
a/Fe classification, and the source of our spectroscopic 
data for oxygen abundance analysis, are given in Table[TJ 

4. DISCUSSION 

4.1. Oxygen in Low- and High-a Halo Stars 

NS10 convincingly showed that halo stars in the so- 
lar neighborhood divide into two groups based on their 
a-element abundances. Their main result is reproduced 
here in FigureQJ,, but note that we show only the 67 stars 
for which we have also derived oxygen abundances. The 
NS10 work is based on the analysis of 94 stars. Neverthe- 
less, our sub-sample is large enough to show the differ- 
ences and similarities between the three groups plotted in 
Figure [TJi. NS10 sorted their stars into high-a halo stars 
(open circles in Figure [I}, low-a halo stars (filled circles), 
and thick-disk members (crosses). The latter were in- 
tentionally included in the work by NS10 for comparison 
with the old disk population, and they were disentangled 
from the halo group by their total Galactic space veloc- 
ity, with thick-disk stars being slower than 175kms . 
While the high-a and thick-disk groups appear to have 
indistinguishable [a/Fe] versus [Fe/H] trends, the low-a 
group, as their name suggests, have systematically lower 
[a/Fe] at any given [Fe/H]. Note, however, that the low 
and high-a groups appear to merge at the lowest [Fe/H] 
values covered by this sample. Although oxygen is an 
a element, it was not included in the NS10 work; their 
[a/Fe] values are based on Mg, Si, Ca, and Ti abun- 
dances. 

The LTE oxygen abundance patterns we derive for the 
NS10 stars are plotted in Figure [TJd. The exact same 
behavior is noted for the three groups, i.e., low-a stars 
have also low oxygen abundances. Note that the [O/Fe] 
abundance ratios appear to increase with lower [Fe/H] 
in a relatively rapid manner, similar to the [a/Fe] case. 
However, this trend is mostly due to non-LTE effects, 
which become important at low [Fe/H] and the some- 
what warmer T e g values of the more metal-poor stars 
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Figure 1. a) a-element abundances as a function of [Fe/H] for the 
thick-disk (crosses), high-a (open circles), and low-a (filled circles) 
halo stars from NS10 included in this work. The two low-a stars 
with large open circles surrounding them are G53-41 ([Fe/H] = 
-1.20) and G150-40 ([Fe/H] = -0.81). b) LTE oxygen abundances 
as a function of [Fe/H] for the stars plotted in panel a), c) Non- 
LTE corrected oxygen abundances. The a-element abundances are 
those derived by NS10; the oxygen abundances are those derived 
in this work (TW). Typical error bars are shown at the bottom 
right side of each panel. 

in this sample (cf. Section [3]). As shown in Figure [It, 
the [O/Fe] versus [Fe/H] relation, corrected for non-LTE 
effects, is nearly flat up to [Fe/H] ~ —0.7 for both the 
low- and high-a groups, and also for the thick disk. A 
hint of a knee towards lower [O/Fe] abundance ratios 
at [Fe/H] ~ —0.7 for the high-a group is detected, and 
it will be confirmed later in this paper when we intro- 
duce additional (literature) data. Hereafter, the oxygen 
abundances used in our discussion are those corrected for 
non-LTE effects, i.e., those shown in Figure [TJ;. 

The star-to-star scatter in [O/Fe] at a given [Fe/H] ap- 
pears to be the largest for the low-a group. Simple linear 
fits to the [O/Fe] versus [Fe/H] relation result in a la 
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scatter of 0.056 dex for the high-a stars, 0.047 dex for the 
thick-disk stars, and 0.075 dex for the low- a group. We 
note, however, that there are two stars which show very 
low oxygen abundances ([0/Fe]< 0.2) and they increase 
significantly the scatter of the low-a group. Excluding 
these stars, the 1 a scatter of the low-a group (0.058 dex) 
is essentially the same as that of the high-a halo stars. 

As pointed out by NS10, and as mentioned before, 
the low- and high-a populations appear to merge below 
[Fc/H] ~ —1.4, making it more difficult to disentangle 
them using chemical abundances. Note, however, that 
the separation is more clear if [a/Fe] is used instead of, 
for example, [Mg/Fe]. Moreover, if we calculate a new 
a element abundance as [a/Fe]' = (4 [a/Fe] + [0/Fe])/5, 
the separation appears even clearer. 

Our sample includes five low-a and two high-a halo 
stars with [Fe/H] < —1.4. With the exception of one low- 
a star, which has an [O/Fe] abundance ratio similar to 
the average of that for the high-a stars at [Fe/H] < — 1.4 
(HD 219617, [Fe/H] = -1.45, [O/Fe] = 0.50), Figure [Tfc 
strengthens the classification suggested by NS10. Al- 
though four of the five low-a halo stars with [Fe/H] < 
— 1.4 have [O/Fe] abundance ratios that are marginally 
consistent, considering the 1-er error, with the average 
[O/Fe] of the most metal-poor high-a stars, it is highly 
unlikely that all, simultaneously, have underestimated 
oxygen abundances. Indeed, the probability that all 
these four stars have [O/Fe] abundance ratios greater 
by 1 a is only about 1 %. Similarly, the probability that 
the two high-a stars both have [O /Fe] lower by 1 a is 
about 10%. Thus, the ambiguity regarding population 
membership at [Fe/H] < — 1.4 may affect only one star: 
HD 219617. If instead of considering this object as a low- 
a star we assume that it belongs to the high-a group, our 
conclusions remain unaltered. For example, the 1 a scat- 
ter values mentioned in the previous paragraph do not 
change by more than 0.001 dex. Therefore, we conclude 
that the "merging" of low- and high-a populations at 
[Fe/H] < —1.4 do not affect our conclusions in a signifi- 
cant manner. 

The two stars with the lowest oxygen abundances in 
Figure [It are G53-41 ([Fe/H] = -1.20, [O/Fe] = 0.18) 
and G150-40 ([Fe/H] = -0.81, [O/Fc] = 0.19). Inter- 
estingly, these stars have perfectly normal low-a abun- 
dances, i.e., although their [a/Fe] abundance ratios are 
low compared to the high-a and thick-disk stars, they 
are not significantly lower than those of a typical low- 
a halo star. NS10 noticed that these stars also have 
very high Na abundances, suggesting that the gas which 
formed them was polluted by nucleosynthesis products 
from nearby asymptotic giant branch (AGB) stars, sim- 
ilarly to Na-enhanced stars in globular clusters (see Sec- 
tion |L4] for more details). 

Not surprisingly, the oxygen and a-element abun- 
dances correlate well, as shown in Figure O with the 
exception of the two stars with the lowest [O/Fe] abun- 
dance ratios mentioned above. Although the star-to-star 
scatter of the [a/Fe] versus [O/Fe] relation is low for 
the high-a and thick-disk stars (1 a scatter of 0.033 and 
0.017 dex, respectively, for a simple linear fit), that for 
the low-a stars is clearly larger (0.057 dex), even if we 
exclude the two stars with [0/Fe]< 0.2 (0.055 dex). 




Figure 2. [a/Fe] versus [O/Fe] relation for the stars in Figure[T] 
Typical error bars are shown at the bottom right corner. 

4.2. Galactic Chemical Evolution of Oxygen: 777 ' nm 
Ol Triplet Analyses of Solar Neighborhood Stars 

In order to put the oxygen abundance trends of low- 
and high-a halo stars, as inferred from the 777 nm Ol 
triplet lines, into a broader context, in Figure [3] we show 
also the non-LTE corrected [O/Fe] abundance ratios by 
M06 and R07. Stars from the NS10 work included in 
these previous studies were excluded; i.e., we used the 
results obtained in this paper instead of the literature val- 
ues. The use of M06 and R07 data allows us to study the 
GCE of oxygen from [Fe/H] = -2.3 to +0.3, and there- 
fore that of the different stellar populations it includes, 
albeit only their solar neighborhood members. We note 
that, although not identical, the stellar parameter de- 
termination and non-LTE corrections applied in M06, 
R07, NS10, and this work, are similar. This prevents 
systematic offsets from biasing our data and artificially 
introducing noise to the chemical abundance trends. As 
shown in Figure |3l stars of common populations but from 
different data sets connect nicely, suggesting that system- 
atic differences between these three studies are indeed 
not very important. 

The stars from M06 are shown with open squares 
in Figure [3] they are all metal-poor main-sequence 
and turn-off s t ars w ith h alo kinemati c s sele cted by 
lAkerman et~aLl (|2004f l and iNissen et ail (|2004f ). These 
authors mention that their sample has halo kinematics, 
without given further details. M06 verified that they 
indeed have halo kinematics based on their large total 
Galactic space velocities (Vtot)- In fact, the Galactic 
space velocities of this sample fulfill the same criteria 
used by NS10 for selecting halo stars, meaning that they 
all have Vtot larger than 180 kms -1 , except one star, 
BD +083095, which has halo kinematics according to its 
admittedly uncertain Hipparcos parallax, but it may ac- 
tually be a thick-disk star based on its spectroscopic par- 
allax. Nevertheless, the uncertain membership of this 
single star does not affect our results. 

Most of the objects from M06 appear to be the natural 
extension of the high-a, high-oxygen abundance popula- 
tion down to [Fe/H] = —2.3. Four objects from M06, 
however, seem to have low oxygen abundance and are 
more likely associated with the low-a population. These 
stars are CD -42 14278 ([Fe/H] = -1.86, [O/Fe] = 0.46), 



6 



RAMIREZ ET AL. 




-2.0 



-1.5 



-1.0 
[Fe/H] 



-0.5 



0.0 



Figure 3. Non-LTE oxygen abundances inferred from the 777 nm Ol triplet as a function of [Fe/H]. Crosses, large open circles, and filled 
circles correspond to the stars plotted in Figure [T] Small open circles and asterisks represent thin- and thick-disk stars, respectively, from 
R07. The membership criterion for the disk stars is based on the abundances, not the kinematics. Hence the perfect separation of thin- and 
thick-disk stars. Open squares and star symbols correspond to halo stars from M06 and R07, respectively. Typical error bars are shown at 
the bottom left corner. 



G24-3 ([Fe/H] = -1.47, [O/Fe] = 0.31), HD 146296 
([Fe/H] = -0.72, [O/Fe] = 0.31), and HD 160617 
([Fe/H] = -1.70, [O/Fe] = 0.30). They can be identi- 
fied in the color version of Figure [3] as red open squares. 

The R07 data includes thin-disk, thick-disk, and a 
few halo stars. Along with the M06 and NS10 stars, 
objects from R07 are all in the solar neighborhood, 
i.e., within a volume with radius of a few hundreds of 
parsecs. This region is dominated by thin-disk stars. 
Studying thick-disk or halo stars in-situ would require 
observing objects at distances of order 1 kpc in the 
direction perpendicular to the Galactic plane, as the 
thick-disk scale heig ht has been claimed to be between 
0.6 an d 1.5 kpc (e.g.. iGilmore fe Reidl[T98l iSieeel et all 
20021: iCabrera-Lavers et al l 120051 : Ide Jong et all 120101: 
Mateu et al. J|201lO. Although a kinematic criterion could 
be applied to disentangle the disk populations in the 
solar neighborhood, it has been shown that there is 
not a perfect one-to-one correspondence with the abun- 
dances (see also below). Thus, for local stars within 
a few hundreds of parsecs of distance, a chemical tag- 
ging approach would be more appropriate. In this 
work, we separate thin-disk and thick-disk stars by 
their [O/Fe] abundance ratios instead of their kinemat- 
ics. A broken line with nodes at ([Fe/H], [O/Fe]) = 
(-0.03, 0.22), (-0.55, 0.35), (-1.00, 0.35) is (somewhat 
arbitrarily) used as the thin/thick disk boundary. In Fig- 
ure [3l our chemically-tagged thin-disk stars are shown 
with small open circles, while the thick-disk members 
are plotted with bold asterisks. The latter also appear 
to be a natural extension of the NS10 data for thick- 
disk stars up to [Fe/H] ~ —0.1. The apparent lack of 
disk stars with [O /Fe] abundance ratios intermediate be- 
tween those of a typical thin-disk and thick-disk star at 
[Fe/H] ~ —0.5 is most likely due to sample selection bi- 
ases (see Section |4~3"|) . 

A few objects with halo kinematics from R07 are shown 
in Figure [3] with five-pointed stars. These stars have a 
probability greater than 50 % of being halo members ac- 
cording to the kinematic criterion employed by R07 (their 
Sect. 3.3). In summary, the thin disk, thick disk, and 
halo populations are assumed to have Gaussian Galactic 



space velocity distributions, wit h mean U, V, W values 
and velocity dispersions gi ven bv ISoubiran et all (|200l 
for the thin/thick disk and lChiba fe Beersl (|2000l ) for the 
halo. 

Only one of the halo stars from R07 seems to be- 
long to the low-a population: HIP 4544 ([Fe/H] = 
-0.81, [O/Fe] = 0.30). We note that, according to 
R07, HIP 4544 has a high probability of being a thick- 
disk member (43%), but its low oxygen abundance 
clearly suggests th at it is a low-a halo object instead. 
iReddv et all (|2006l ) have measured the [a/Fe] abundance 
ratio for this object, which is 0.23, a value that at 
[Fe/H] = —0.81 appears too high for a low-a halo star, 
but it is only marginally consistent with the thick-disk 
trend, so it is not at all clear w hat q-element popul ation 
this star belongs to. Note that IReddv et al.l (|2006l ) and 
NS10 a-element abundances are not necessarily on the 
same scale, so systematic differences between these two 
studies could be responsible for this apparent discrep- 
ancy. Combining all three data sets, we find that the 
high-a population shows [O/Fe] abundance ratios that 
decrease slightly from [O/Fe] ~ 0.60 at [Fe/H] = -2.3 to 
O/Fe] ~ 0.55 at —0.7, and from there more abruptly to 
O/Fe] ~ 0.35 at [Fe/H] ~ -0.1. 

The chemically-tagged stellar populations from Fig- 
ure [3] have been re-plotted in Figure 31 in three panels: 
a) disk, b) high-a, high-oxygen halo, and c) low-a, low 
oxygen halo. In panel a), the solid line which has a knee 
at [Fe/H] ~ —0.6 is a fit to the high-a disk data (the 
"thick-disk" ) whereas the solid curve is a cubic fit to the 
low-a disk data (the "thin-disk"). These fits are repro- 
duced in panels b) and c). In panel b), the dashed line 
which has a knee at [Fe/H] ~ —0.7 is a fit to the high-a 
halo data, and is reproduced in panel c). The "broken" 
linear fits were made by allowing the knee location to 
be a free parameter. We monitored the 1 a star-to-star 
scatter around these linear fits to make sure that these 
relations closely correspond to a minimum standard de- 
viation. We also ensured that the fits resulted in two 
linear segments connecting smoothly at the knee. 

The cubic fit to the thin-disk data has a star-to-star 
[O/Fe] scatter of only 0.04 dex, which is compatible with 
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Figure 4. a) Oxygen abundance trends of thin- (circles) and 
thick-disk (crosses and asterisks) stars. The solid curve and broken 
line are fits to the thin- and thick-disk data, respectively, and they 
are shown also in panels b) and c). b) Oxygen abundance trend 
of high-a halo stars. The dashed broken line is a fit to the high-a 
halo data and it is shown also in panel c). c) Oxygen abundance 
pattern of low-a halo stars. Typical error bars are shown at the 
bottom left corner of each panel. 

a zero cosmic scatter, i.e., it suggests that observational 
errors alone explain it. This probably reflects the fact 
that the thin-disk stars have been born from well-mixed 
material at late stages in the history of Galactic evolu- 
tion. The steep decline in [O/Fe] abundance ratio with 
increasing [Fe/H] has often been attributed to the chem- 
ical pollution of the interstellar medium by Type la su- 
pernovae (SNIa), which dominate GCE only after a few 
billion years since the birth of the Galaxy. 

The high-a halo and thick-disk populations appear 
very similar, but the linear fits suggest a small down- 
wards offset for the thick disk relative to the halo. Also, 
the location of the knee seems to be about 0.1 dex lower 
for the high-a halo group. The star-to-star scatter of 



these fits is about 0.06 dex. Although observational er- 
rors are larger for these groups of stars, mainly because 
they are more distant and their atmospheric parameters 
cannot be determined as precisely as those of most thin- 
disk stars, some of this scatter could be real. Estimates 
of [O/Fe] errors are about 0.05 dex for most of these ob- 
jects. We note also that the star-to-star scatter of the 
thick-disk stars with [Fe/H] < —0.6 is about 0.07 dex, 
suggesting that a less well-mixed gas, or more likely a 
mixture of different gases, gave origin to these objects 
compared to the high-a halo stars. 

The fact that the high-a halo and thick-disk star 
oxygen abundance trends are best fit with a broken 
line suggests that their more metal-poor members were 
born from gas enriched mainly by the pollution of 
massive stars, whose yields (at the end of their lives 
as Type II supernovae, SNII) have large O/Fe abun- 
dance ratios. The shallow negative slope of the [O/Fe]- 
[Fe/H] relation at low [Fe/H] (< -0.65) could be fully 
explained by metallicity-dependent SNII yields (R07). 
Later ([Fe/H] > —0.6), SNIa contributed significantly 
to the chemical evolution of these populations, quickly 
lowering their [O/Fe] ratios while [Fe/H] continued to 
increase. That the high-a halo stars have a knee at a 
slightly lower [Fe/H] is a possible indication that the 
transition happened earlier in the halo than in the thick 
disk, or that the star formation rate was somewhat slower 
in the halo. Note that the [Fe/H] errors are about 
0.05 dex in this compilation of oxygen abundances, while 
the difference in [Fe/H] for the location of the high-a 
halo and thick-disk knees is about 0.1 dex, implying that 
this difference, although marginally consistent with zero 
within the errors, is unlikely due to observational uncer- 
tainties. 

As shown in Figure Ht, the low-a halo stars do not 
follow any of the other [O/Fe] trends. They appear as 
a continuation of the thin disk, but this connection is 
very weak on light of other well-known observational ev- 
idence such as the stellar kinematics and age (thin-disk 
stars are younger and have cold kinematics, the complete 
opposite of what is normal in halo stars). On average, 
and excluding the stars with very low oxygen abundances 
discussed above, these stars have [O/Fe] abundance ra- 
tios that are about 0.2 dex lower than those of high-a 
and thick-disk stars. There seems to be a downwards 
trend of [O/Fe] with increasing [Fe/H], but the star-to- 
star scatter is large. In any case, the mean [O/Fe] in 
the low-a halo population decreases from [O/Fe] ~ 0.5 
at [Fe/H] ~ -1.8 to [O/Fe] ~ 0.35 at [Fe/H] ~ -0.7. 

Low [O/Fe] abundance ratios at low [Fe/H] have been 
typically attributed to populations with a history of slow 
star formation rate. They have been observed in the 
dwarf satellite galaxies of the Milky Way, suggesting 
that the low-a, low-oxygen halo stars are remnants of 
merger processes that occurred early in the history of 
our Galaxy, or perhaps stars that have been stripped 
from their parent satellite galaxies as they came close or 
passed through the solar neighborhood in their Galactic 
orbits. 

4.3. Kinematics and Chemical Abundances 

A very important tool to explore further the GCE in- 
terpretations is the Toomre diagram, which is shown for 
our sample of stars in Figure[5l Here we plot the Galactic 



RAMIREZ ET AL. 



400 



300 




100 







-500 -400 -300 -200 
V (km s" 1 ) 

Figure 5. Toomre diagram for the stars plotted in Figure [3] 
Dashed lines correspond to speeds of 65 and 180 km s -1 . 

space velocity V (which goes in the direction of disk rota- 
tion and is measured relative to the Sun), and the other 
two components combined: UW — VU 2 + W 2 (U is de- 
fined as the Galactic space velocity component towards 
the Galactic center while W is the component perpen- 
dicular to the plane). This diagram is sometimes used 
to separate stars kinematically into thin-disk, thick-disk, 
and halo members. Although such kinematic member- 
ship criteria have been proven useful, particularly for 
solar neighborhood stars, as larger samples of stars are 
analyzed using better spectroscopic data and performing 
more careful elemental abundance work, we are starting 
to see that kinematics and chemical composition do not 
make a one-to-one correspondence for the local stellar 
populations. For example, although thin-disk stars tend 
to have cooler kinematics (i.e., V closer to zero and low 
UW) compared to thick-disk members (lagging V ve- 
locities and higher UW), there is an important number 
of stars with chemical composition typical of that of a 
thin-disk (thick-disk) star which have thick-disk (thin- 
disk) kinematics. The observed fraction of these stars 
with "ambiguous" kinematics and chemical abundances 
cannot be fully accounted for by observational errors 
(Ramirez, Allende Prieto, & Lambert, in preparation). 

Figure [5] shows that the chemically-tagged thin-disk 
stars, as a group, tend to have cold kinematics, i.e., they 
rotate fast (V closer to zero) and do not depart much 
(small UW) from the mid-plane of the Milky Way's disk. 
Thick-disk stars, on the other hand, rotate slower (i.e., 
lag the Sun and the thin-disk stars, as a group) and have 
larger U an d W velocities. This has been known f or many 
years (e.g.. lSoubirarilT99 3; Sou biran et alJl2003l h and it 
has been used to select one group o f stars or the other 
in chemical abund ance studies (e.g.. iBensbv et aDl2004t 
iReddv et al.ll2006h . However, there is significant overlap 
between the two groups, and because of that, kinematic 
selection of thin/thick disk stars has often avoided the in- 
termediate region, leading to thin/thick disk abundance 
patterns that may be affected by severe kinematic bi- 
ases. An extended discussion of this possible bias will 
be given in Ramirez, Allende Prieto, and Lambert (in 
preparation). Here, we note that, although the thin- 
disk and thick-disk groups appear on average to separate 
kinematically, there is a significant number of stars with 



thin-disk abundances but thick-disk kinematics, and vice 
versa. If we make a simple kinematic separation, with a 
boundary of total speed equal to 65kms _1 , which corre- 
sponds to the inner dashed line shown in Figure[5j we find 
that these stars with ambiguous kinematics and elemen- 
tal abundances amount to about 30 % of each sample, 
a number that is certainly not negligible, as has been 
sometimes assumed or ignored altogether. 

The thick disk has long been thought to be formed 
from a violent me rger event early i n the history of 
the Galaxy (e.g., iQuinn et afl [1993). This merger 
would have destroyed a previously formed disk, heat- 
ing the stars into more eccentric orbits. It is sug- 
gested by some authors that thick-disk stars could have 
been formed in the merging galaxies a s well as in 
the original disk (e.g.. [Abadi et al.1 120031 : iBrook et al.1 
120051 iKobavashi fc Nakasatd 120111 ). These scenarios, 
however, are not able to fully explain the large frac- 
tion of stars with ambiguous kinematics and abun- 
dances. An alterna t ive pi cture, revived by the work of 
iSellwood fe; Binnevl (|2002[ ). involves internal processes, 
particularly radial mixing, which ex plains many of the so- 
lar neighborhood observa bles ( e.g., iSchonrich fc Binnevl 
I2009al|bl : lLoebman et~aT1 120111 ). Nevertheless, recent 
SEGUE/SDSS observations of the Galaxy on a larger 
scale appear to be incons istent with th ese models 
(e.g., Schlcsin ger et al.l 120121 see, however, IBovv et all 
12012a HQ). The issues of sample completeness and sam- 
ple selection functions need to be fully addressed be- 
fore observations of solar neighborhood stars are used 
to determine which scenario is more realistic. Up- 
coming and ongoing hig h-resolution sp ectroscopic sur- 
yeys such as HERMES (iFreemanl |2010|) and AP OGEE 
(|Allende Prieto et al.ll2008l : lMaiewski et . al.ll2010f ) should 
allow us to solve this problem within the next decade. 

Regarding the halo, there is not a one-to-one corre- 
spondence in chemical abundances and kinematics of 
high-a and low-a stars either, but an average trend was 
detected by NS10. They found that the low-a stars tend 
to have very low V velocities, in fact most of them are in 
retrograde Galactic orbits. Four of the five low-a stars 
we found in M06 and R07 and which are not included 
in NS10 also have retrograde orbits, strengthening their 
conclusion. In their update of the R07 work, Ramirez, 
Allende Prieto, and Lambert (in preparation) also find 
a number of additional halo stars in retrograde orbits 
which have low oxygen abundances relative to the "nor- 
mal" halo. 

In Figure [51 our oxygen abundances are plotted against 
the orbital parameters r max (maximum distance from the 
Galactic center), |2:| max (maximum height with respect to 
the Galacti c plane), and e-ma* (ma ximum eccentricity) , as 



derived by iSchuster et al.l (|2012l ) from 5 Gyr orbit inte- 
grations computed using the NS10 star's kinematic data 
and a detailed semi-analytic model for the Milky Way 
potential. The orbital paramete rs plotted in F i gure E] 
correspond to those obtained by ISchuster et al.l ()2012f ) 
with a realistic non-symmetrical Galaxy model. Note 
that the stars plotted in Figure [6] are those from Fig- 
ure Q] and that the symbols used there are the same as 
those used in Figure [T] as well. 

All of the high-a, high-oxygen halo stars have orbits 
that do not go beyond 20 kpc of distance from the Galac- 
tic center, contrary to about half of the low-a, low- 



OXYGEN IN LOW- AND HIGH-ALPHA HALO STARS 



9 



40 



o. 30 



1 20 



15 

10 

5 


1.0 

0.8 

0.6 

0.4 

0.2 



» 



+ 



m 



® 



o 



°o o o 

° Oo CP 



® 
® 



-H 

O 

O O 

" P o 



oo 



o° O 



0.2 



0.3 



° o 




+ + + + 


+ 


+ + 


+ 


+ 

. . . 1 . . . 


+ 


0.5 


0.6 



0.4 
[0/Fe] 

Figure 6. Orbital parameters r max (maximum horizontal dis- 
tance from the Galactic center), z max (maximum vertical distance 
from the Galactic plane), and e max (maximum eccentricity) as a 
function of [O/Fe] for the NS10 stars analyzed in this work. Sym- 
bols are the same as in Figure [l] 

oxygen stars, which have r max up to about 40kpc. Sim- 
ilarly, the high-a, high-oxygen stars reach heights above 
the Galactic plane up to about 8 kpc whereas some low- 
a stars have orb its that go about tw ice as high (see 
also Figure 8 in ISchuster et al.l 120121) . The orbits of 
thick-disk members, on the other hand, are smaller in 
both directions, and they are also less eccentric than 
those of both types of halo stars. In fact, a major- 
ity of the latter have e max greater than 0.8, whereas 
thick-disk stars have e max ~ 0.35 ± 0.15. Note also 
that there are more high-a, high-oxygen stars with e max 
lower than 0.8 than low-a, low-oxygen stars. The im- 
plications of these distinct orbital distribu tions were al- 
ready discussed by ISchuster et al.l (|2012l ). In particu- 
lar, they reinforce the idea that the two halo popula- 
tions require different formation scenarios, with the low- 
a group being accreted stars. Interestingly, the distribu- 
tion of r max values for globular clusters is fully consistent 
with the scatter of r max values se en in Figure El for the 
low- en, low-oxygen h alo stars (e.g.. iDauphole et al.lll996t 
iDinescu et al 



ygen h-c 
. 1999). 



4.4. Field Halo Stars Born in Globular Clusters 

Based on their kinematics, NS10 proposed that the 
low-a halo stars could have been born in the dwarf satel- 



lite galaxies of the Milky Way, with some of them prob- 
ably originating from the globular cluster (GC) wCen. 
The more detailed chemical composition ana lysis made 
by these authors in iNissen fc Schusterl (|2011l ). however, 
revealed more differences (e.g., in a, Na, and Ba/Y) than 
similarities (e.g., in Ni and Cu) between the low-a halo 
stars and uj Cen. As suggested by them, perhaps chemi- 
cal evolution in uj Cen was different for its inner (or more 
bound) and outer (less bound) regions, explaining the 
present-day differences. 

The association of groups of field halo stars with uj Cen 
based on chemical analysis is tempting, as the many ex- 
amples that can be found in the literature demonstrate, 
including the NS10 work. We must be reminded, how- 
ever, that of all of the Milky Way's GCs, uj Cen is the 
most comple x example, exhibiting a wide range of stel- 
lar ag es (e.g. . lHugh.es fe Wallerst cin 2000; [Stanford et~al 
2006) and metallicities (e.g., iNorris fc Da Costal 11995 



Frinc habov et al.l 12002). The latter imply that chemi 



cal evolution within the cluster has occurred following 
not a single but a number of episodes of star forma- 
tion. Indeed, large chemical abundance surveys of uj Cen 
stars suggest distinct chemical evolution paths followed 
by a number of clearly identified sub-populations (e.g., 
Uohnson fc Pilachowskill201tt iMarino et al.ll2011b[ ). 

One of the most notable chemical properties of 
u> Cen, observed also in most oth er GCs, is the so- 
called Na-0 anti-correlatio n (e.g . . INorris fc Da Costal 
|1995bt iGratton et al.1 12001L 120071: iCarretta et all 1 2009 1 : 
Johnson fc Pilachowskil l2010t iD'Antona et al.l boill : 
iMarino et al.l 1201 lbl ). a property that is not seen in 
field halo stars. It has been suggested that this anti- 
correlation is due to in-situ mixing of inte rmediate-mass 
AGB star nucleosynthesis products (e.g., iVentura et al.l 
l2001t iGratton et al.ll2004f ) . With our oxygen abundance 
data and the Na abundances from NS10 we can now ex- 
plore another possible connection in the form of the Na-0 
anti-correlation. 

In Figure[7]we plot [Na/Fe] versus [O/Fe] for the stars 
studied in this work. With the exception of G53-41 and 
G150-40, the low-a, low-oxygen abundance halo stars all 
have sub-solar [Na/Fe] abundance ratios. This observa- 
tion is consistent with the o xygen and sodium abundanc e 
data for uj Cen red giants by INorris" fc Da Costal (11995bD . 
A similar conclusion can be r eac hed by lookin g at the 
Uohnson fc Pilachowskil (|2010f l or IMarino et all l|2011bft 
larger data sets. We note, however, that the spread of 
oxygen and sodium abundances measured in u) Cen stars 
is large, and in fact the [Na/Fe] versus [O/Fe] relation 
in uj Cen overlaps also with the location of the high-a 
halo and thick-disk stars studied in this work. Only the 
presence of G53-41 and G150-40, the stars with lowest 
[O/Fe] and highest [Na/Fe] in our sample, hints at an 
Na-0 anti-correlation for the low-a group. In fact, their 
peculiar chemical composition can be attributed to pol- 
lution by nearby AGB stars to the proto-stellar gas, in a 
similar fashion to globular cluster stars (NS10), although 
note that the abundance anomali es could also be due t o 
fast rotating massive stars (e.g., Decrcssin et al. 2007). 
In any case, it is clear that the majority of low-a stars 
do not exhibit an obvious Na-0 anti-correlation, which 
in principle further weakens the uj Cen connection. Nev- 
ertheless, it is important to point out that there are no 
stars with this type of peculiar composition in the high-a 
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Figure 7. [Na/Fe] versus [O/ Fe] relation for the s tars in Figure[T] 
Sodium abundances are from Nisscn & Schuster (2010). Typical 
error bars are shown at the bottom left corner. 

halo or thick-disk groups. 

The stars G53-41 and G150-40 are the only field stars 
showing the classical signatures of abundance anomalies 
in second generation GC stars, i.e., enhanced Na and 
depleted O. Other abundance peculiarities are also dis- 
cussed in Section l4~5l Since we have analyzed 67 stars, we 
conclude that the fraction of metal-poor field stars orig- 
inating from second-generation GC stars is about 3 %. 
Adopting a binomial distribution, which is appropriate 
in this case given the relatively low number of objects 
and the fact that there are two possible "outcomes" for 
each star, i.e., field and GC, an error bar can be esti- 
mate d from the varia nce of the probability distribution 
(e.g.. lBevTn gton 1969, Chapter 3): a 2 — np(l — p), where 
n = 67 is the number of stars and p the probability of 
"success" (p = 2/67 = 0.03). We find a = 1.4, which 
implies a probability error of 1.4/67 = 2%. 

Of course, the actual fraction of halo field stars orig- 
inally formed in GCs may be significantly higher than 
the value of 3 ± 2 % derived above, as the clusters have 
likely also contributed to the halo field with "normal" 
(i.e., first generation) stars, which may be hard to dis- 
tinguish from the bulk of halo field stars observed today. 
Thus, our oxygen abundances and the Na abundances 
from NS10 suggest that the fraction of halo stars born in 
globular clusters is at least 3±2 %0 Indeed, although the 
fraction of field metal-poor giants with anomalous CN 
and CH bands (typica l of second-generation GC stars) is 
only 3% according to iMartell etaTJ (|201lD . a minimum 
of 17% of the present-day mass in the halo field origi- 
nated from GCs. Us ing the binomial distribution on the 
IMartell et all (|201l! ) results, we obtain a more precise 
value of 2.85 ± 0.70% for the first of these fractions, in 
excellent agreement with our estimate of 3 ± 2 %. Note, 
however, that our highly reliable O and Na abundances 
connect these objects to GCs in a more direct way than 
the intensities of CN and CH bands. 

A comparison between the two peculiar low-a stars 
(G53-41 and G150-40) and the Na and O abundances 
observed in second generation stars in the globular clus- 
ter M71 is instructive, as this cluster has [Fe/H] = —0.8 
(|Melendez fc Cohen! [20091 ). In relation to first gcncra- 

9 Errors in our [O/Fe] and NSlO's [Na/Fe] abundance ratios are 
too small to have a significant impact on this lower limit. 



tion stars, the second generation stars in M71 have O 
depleted by on ly about 0.1 dex while Na is enhanced by 
about 0.4 dex (iMelendez fc Cohenl [20091 . Similarly, for 
M4 ([Fe/H] = -l.D. lMarino et al J ([20111 find an O de- 
pletion of about 0.25 dex and a Na enhancement of about 
0.4 dex. These abundance variations are fully compatible 
with the ones observed between G53-41 and G150-40 and 
the bulk of low- a halo stars (Figure [7|) . While oxygen in 
these two dwarfs is depleted by about 0.2±0.1 dex, Na 
seems enhanced by about 0.4±0.1 dex. This quantitative 
comparison reinforces the idea that the two anomalous 
dwarfs may have been formed in GCs. If that is the case, 
G53-41 and G150-40 would be the first field stars with 
firm O/Na signatures of being originated in GCs. The 
discovery of these signatures certainly shows the benefit 
of high precision differential abundance studies. 

4.5. Oxygen Compared to Other Elements 

It is interesting to note that, even though the high-a, 
low-a, and thick-disk stars studied in this work sepa- 
rate in both [Na/Fe] (Figure 6 in NS10) and [O/Fe] (our 
Figure [1} when each of these abundance ratios is plot- 
ted against [Fe/H], the run of [O/Na] with [Fe/H] is es- 
sentially indistinguishable between the three groups, as 
shown by Figure [5] The obvious exceptions are, again, 
G53-41 and G150-40. This striking similarity in elemen- 
tal abundance ratios [O/X] between the three groups is 
seen also for the a-elements as well as Zn. Other elements 
show systematic offsets in the distribution of [O /X] abun- 
dance ratios of the three stellar populations, although for 
Y and Ba the star-to-star scatter is too large to distin- 
guish them, if present. 

In Figure the range in the [O/X] axis has been set 
equal to 1.0 dex for all panels, except that for Ba, for 
which it is 1.2 dex. This allows a fair comparison of the 
star-to-star scatter differences and offsets between high- 
a halo stars, low-a halo stars, and thick-disk stars. For 
example, it is clear that the run of [O/Zn] with [Fe/H] is 
the tightest of all, especially if we exclude the three low-a 
stars with the lowest [O/Zn] . Two of these three stars are 
G53-41 and G150-40, which have been mentioned before 
as having peculiar composition. The other object, which 
is the one with the lowes t [O/Zn] in our sample, is Gil 2- 
43. As pointed out by iNissen fc Schuster! (|2011f) . this 
object is one of the components of a wide binary, and, 
along with its companion (with which it has consistent 
chemical abundances), they seem to also have peculiar 
abundances (in [Mn/Fe] and [Zn/Fe]) when compared to 
the mean low-a trends. 

Excluding G53-41, G150-40, and G112-43 from Fig- 
ure [8] results in tighter correlations and stronger pop- 
ulation similarities for [O/a], [O/Na], and [O/Zn]. In 
addition to the latter abundance ratios, the outlier stars 
appear to be very peculiar in [O/Y] and [O/Ba]. The 
most extreme case is [O/Na], but, interestingly, only for 
G53-41 and G150-40, i.e., not for the binary component 
G112-43. For [O/Y] and [O/Ba], only one of the two 
stars with very low oxygen abundance (G53-41) has an 
unusually high (relative to other low-a stars) [Ba/Fe] 
abundance ratio ([Ba/Fe] = +0.24 according to NS10), 
whereas the other (G150-40) shows marginally low [O/Y] 
and [O /Ba] which may be considered normal for a low-a 
halo star. Barium abundance anomalies relative to nor- 
mal field halo stars have also been seen in u Cen members 
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(e.g.. INorris fc Da Costal [19951 ISmith. et all 120001) . but 

the enhancement observed in those stars is significantly 
higher than that measured in G53-41. 

5. CONCLUSIONS 

Non-LTE oxygen abundances from the 777 nm O I 
triplet lines have been derived for as many as possible of 
the stars in the work by NS10. These authors have de- 
rived very precise atmospheric parameters and elemental 
abundances (excluding oxygen) for their sample stars, al- 
lowing them to clearly separate the field halo stars into 
low- a and high- a groups. 

We find the run of [O/Fe] abundance ratios with [Fe/H] 
of high-a halo and thick-disk stars very similar, while 
that of low- a halo stars is systematically lower by about 
0.2 dex and it has, in general, a larger star-to-star scatter 
compared to the other two groups. A few additional 
low-a, low-oxygen abundance halo stars are identified in 
previously published works. Their kinematic properties 
strengthen the hypothesis by NS10 that these objects 
may have originated in dwarf satellite galaxies early in 
the history of the Milky Way. A connection between 
the low-a, low oxygen halo stars and u> Cen is not well 
established, unless assumptions about the early chemical 



stars in Figure [T] Elemental abundances other than that of oxygen are 
1.0 dex in all panels except [O/Ba]. 

abundance distribution within this extremely complex 
globular cluster are made. 

Our oxygen abundance data for the three groups of 
stars studied by NS10 exhibit a behavior that is simi- 
lar to that of the a-elements. The exceptions are two 
stars, G53-41 and G150-40, which seem to be the first 
firm candidates of field halo stars born in globular clus- 
ters, although probably not uj Cen, which has been pre- 
viously argued as one of the main contributors of low-a 
field halo stars. Both G53-41 and G150-40 show the clas- 
sic signatures of abundance anomalies in globular clus- 
ter stars, namely very low oxygen and highly enhanced 
sodium abundances. Since these properties are seen in 2 
of the 67 stars studied in this work, we estimate that the 
contribution of globular clusters to the local field metal- 
poor (—1.6 < [Fe/H] < —0.4) stellar population is at 
least 3 ±2%. 
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Table 1 

Stellar Parameters and Oxygen Abundances. 1 
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1 The atmospheric parameters T c ff,logg, [Fe/H] and o/Fe flags are from Nisscn & Schuster (2010). 



